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Abstract — The article proposes a heuristic approximation ap- 
proach to the bin pacliing problem under multiple objectives. In 
addition to the traditional objective of minimizing the number 
of bins, the heterogeneousness of the elements in each bin is 
minimized, leading to a biobjective formulation of the problem 
with a tradeoff between the number of bins and their heteroge- 
neousness. An extension of the Best-Fit approximation algorithm 
is presented to solve the problem. Experimental investigations 
have been carried out on benchmark instances of different size, 
ranging from 100 to 1000 items. Encouraging results have been 
obtained, showing the applicability of the heuristic approach to 
the described problem. 



I. Introduction 

THE term "bin packing" describes a class of well-known, 
classical problems with numerous applications in logis- 
tics, operations research and related disciplines. From single 
dimensional to multi-dimensional problems, various types can 
be identified in practice. Common to all is the overall task of 
packing a finite number of n items into a minimum number of 
bins (knapsacks) subject to a set of practical constraints and 
requirements. These include given capacities of the bins, but 
also other considerations such as irregularly shaped bins, load 
balancing of the bins, etc. 

Numerous approaches including exact, heuristic, and meta- 
heuristic algorithms have been proposed for the resolution 
of bin packing problems, and a rich literature on packing 
problems exists, with important classifications by Dyckhoff 
[1] and more recently Wascher et al. [2]. While the majority 
of approaches is dedicated to single-objective models, only 
minimizing the number of bins used, the multi-objective nature 
of many of these problems becomes more and more obvious. 
Following early work of WASCHER [3], modern heuristics 
such as Particle Swarm Optimization have recently been ap- 
plied to a multi-objective variant of the two-dimensional bin- 
packing problem [4]. For the here considered multi-objective 
bin-packing problem however, no corresponding studies have 
been carried out to our knowledge. 

The remainder of the article is organized as follows. Sec- 
tion II describes the multi-objective bin packing problem and 
its' underlying practical application. A heuristic approximation 
approach is presented in the following Section III. In brief, 
we propose an extension of the well-known best-fit heuristic, 
allowing the computation of a set of solutions that constitute an 
approximation to the set of efficient solutions. Experimental 
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results of the approach to the problem are reported in Sec- 
tion IV, and conclusions are given in Section V. 

II. Problem description 

We consider a bin packing problem where a given number 
of n items has to be packed into n bins, each of capacity c [5]. 
Each item j is characterized by a weight Wj and an additional 
attribute aj. While the weights refer to the size of the items 
and therefore have to be taken into consideration when filling 
up a bin to at most its' capacity c, the attributes aj describe 
properties of the items on a nominal scale. On the basis of 
this description, a comparison of two items i,j is possible 
such that they are either identical with respect to and aj, 
Ui — aj or not: 7^ aj. The goal of packing the items into 
bins can then be modeled as follows. 
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1 if bin i is used 

otherwise 

1 if item j is assigned to bin i 
otherwise 

Ui counts the number of distinct attributes in bin i 

Expression (1) minimizes the number of bins. The second 
objective given in (2) minimizes the average heterogeneous- 
ness of the bins. To do this, the number of distinct attributes Ui 
is counted for each bin i. Unused bins {yi = 0) have a value 
of Ui = 0. Used bins (yi = 1) have a possible minimum value 
of Ui — 1. This is the case when all items in the particular 
bin have the identical nominal attribute. The values of Ui are 
bounded by either the number of items assigned to a bin or 
the number of distinct attributes over all items i. 



Intuitively, the two objective functions are of conflicting 
nature. While a large number of bins allows the packing of 
bins which are each fully homogeneous, leading to a Z2 = 1, 
a solution being minimal for z\ will require the packing of 
items i,j of different Oj and into the same bin. It can 
therefore be suspected that not a single solution x exists in 
the set of feasible solutions X that equally minimizes both 
objective functions z\ and z-2. In brief, this leads to a vector 
optimization problem in which a solution a; G X is evaluated 
with respect to a vector Z{x) = {z\{x)^ -^2(2;)). The resolution 
of the problem has consequently to be seen in the identification 
of all efficient outcomes or the Pareto set P, introduced in the 
following Definitions 2.1 and 2.2. 

Definition 2.1 (Dominance): A vector Z{x) is said to dom- 
inate Z{x') iff Zk{x) < Zk{x')\/k = l,...,K A3k\zk{x) < 
Zfc(x'). We denote the dominance of Z{x) over Z{x') with 
Z{x) di Z{x'). 

Definition 2.2 (Efficiency, Pareto-optimality): The 
vector Z{x),x £ X is said to be efficient iff 
-y3Z{x'),x' e X I Z{x') d Z{x). The corresponding 
alternative x is called Pareto-optimal, the set of all Pareto- 
optimal alternatives Pareto-set P. 

Numerous practical apphcations of the formal model exist. 
In many cases, a minimum number of bins should be used 
when packing a given set of items, however assuring a max- 
imum possible homogeneousness of the items being packed 
together into a single bin. Applications include the storage of 
goods, the storage of music/video data on optical discs, etc. 

III. A HEURISTIC APPROXIMATION APPROACH 
As already mentioned in Section I, numerous algorithms 
have been proposed to solve the single-objective variant of the 
bin packing problem. One important heuristic is the Best-Fit 
algorithm, important for both for its' time complexity as well 
as for its' worst-case complexity [6]. Best-Fit subsequently 
assigns items to the feasible bin having the smallest residual 
capacity. If no such bin exists, the item is assigned to a 
previously unused (new) bin. 

Unfortunately however, Best-Fit only takes into considera- 
tion the weights Wj of the items and the residual capacities 
of the bins when selecting the 'best-fitting' -bin. In order to 
address the problem described in Section II, a method of 
controlling the heterogeneousness of the bins needs to be 
included in the algorithm. Algorithm 1 describes such an 
attempt, allowing the successive computation of alternatives 
with different heterogeneousness levels and therefore provid- 
ing an idea of how to compute an approximation to the vector 
optimization problem given in Section II. 

The modified Best-Fit algorithm is based in principle on 
the conventional method. However, in order to control the 
heterogeneousness of the bins, an additional control parameter 
Umax is used as described in step 8 of Algorithm 1. Starting 
with an initial value of Umax = u = 1, only Best-Fit-bins are 
allowed which are fully homogeneous. This means that an item 
may only be assigned to a bin containing other elements of 
identical attributes aj. In this stage of the algorithm, solutions 
are computed that lead to the lowest possible value of 22 = 1 
as aU bins contain homogeneous items. 



Algorithm 1 Multi-objective Best-Fit algorithm 
Require: s, m 
1: Compute the maximum possible heterogeneousness of a 
bin u 

2: Set w = 1 

y papprox _ 

4: repeat 

5: for m = 1 to m do 

6: Construct a new solution x: 

7: for all n items do 

8: Compute the maximally allowed heterogeneous- 

ness Umax of the Best-Fit-bin: 
Umax = [u\ with probability 1 — (u — LH ) 
Umax = \u] with probability u — [u\ 
9: Compute the Best-Fit-bin with respect to the ran- 

domly determined Umax 
10: Assign i to the Best-Fit-bin 

11: end for 

12: Update P«pp™^ with x: 

Remove all elements in P'^pp^°^ which are dominated 
by Z(xy, 

Add X to p<^pp^°^ if Z{x) is not dominated by any 

element in P'^pp''"^ 
13: end for 
14: Set U — u-\- s 
15: until u >u 
16: Return P°pp™^ 



With increasing value of u, incremented by s in step 14 of 
the algorithm, Best-Fit-bins become possible that have a higher 
heterogeneousness. This concept is randomized throughout the 
generation of the solutions, allowing a gradual transition from 
Umax = u = 1 to the maximum possible heterogeneousness 
Umax = u. Due to the randomness in the algorithm, different 
runs lead to different outcomes. We therefore propose to 
compute a number of solutions with each setting of u, given 
as control parameter m. 

Throughout the algorithm, an archive P^pp^"'^ of the best 
solutions is kept which is returned after the algorithm termi- 
nates. This archive represents an approximation to the true 
Pareto-set P. 

IV. Experimental investigation 

A. Generation of test instances and experimental setup 

In order to test the effectiveness of the multi-objective 
extension of the Best-Fit approximation algorithm, four multi- 
objective test instances have been computed with values of 
n = 100, n = 200, n = 500, and n = 1000. The data has 
been derived taking ^ bins, each of capacity c = 1000, and 
randomly splitting the capacity into five items j,. . . such 
that the weights Wj of the items add up to c: Y^^j^^ Wj = c. 
This means that the so constructed instances have a solution 
for which the minimum number of bins is and therefore 

equal to the trivial lower bound — — - . 

The items of each instance have been randomly assigned 
nominal attributes from a set of five different attributes. Each 



attribute has been selected with equal probability of 0.2. 

We tested the randomized Best-Fit algorithm from Sec- 
tion III on the proposed benchmark instances with control 
parameters s = 0.1 and m = 100. This means that the 
probabiUty of selecting a Best-Fit-bin of maximum hetero- 
geneousness Umax = \u\ over a Best-Fit-bin of maximum 
heterogeneousness Umax = \u\ increases in each step 14 of 
Algorithm 1 with 10%. The inner loop of the algorithm com- 
putes m = 100 solutions with each setting of the parameters. 

In addition to the multi-objective Best-Fit algorithm, a 
Random-Fit algorithm has been tested for comparison rea- 
sons. This algorithm randomly selects a bin that allows the 
assignment of the currently considered item with respect to 
the chosen maximum heterogeneousness Umax- Apart from 
that aspect, the algorithm is implemented identically to the 
pseudo-code given in Algorithm 1. 

Three different input sequences of the items have been 
tested: 

• Sorting the items in decreasing order of Wj 

• Sorting the items in increasing order of Wj 

• Arranging the items in random order 

B. Results 

The experimental investigations revealed that only few effi- 
cient outcomes exist for the instances. Instead of plotting the 
outcomes in figures, we chose to give the data of all found best 
vectors Z{x) = {zi{:r). Z2{x)). The following Table 1 shows 
the results for the smallest instance with n = 100. It can be 
seen, that both Best-Fit and Random-Fit perform comparably 
good given a decreasing or random order of the items. 

TABLE I 

iDENTinED VECTORS {z\ , 22) FOR THE INSTANCE WITH n = 100. THE 
BEST OUTCOMES ARE HIGHLIGHTED BOLD. 



Item order 


Best-Fit 


Random-Fit 


Decreasing Wj 


(22, 


1.000) 


(22, 1.000) 




(21, 


2.048) 


(21, 1.952) 


Increasing Wj 


(25, 


1.000) 


(25, 1.000) 




(24, 


1.108) 


(24, 1.125) 


Random 


(22, 


1.000) 


(22, 1.000) 




(21, 


1.190) 


(21, 1.952) 



Similar results have been obtained for the instance with 
n = 200 as shown in Table II. Again, Best-Fit and Random- 
Fit lead to the best results given an order of the items with 
decreasing Wj. While both identify the vector (43, 1.000), they 
are incomparable with respect to the other best found outcomes 
(42, 1.214) and (41, 1.902). 

For the next instance with n = 500, given in Table III, 
Random-Fit appears to lead to superior results, however with 
a very small distance to Best-Fit. While the assignment of 
items in increasing order of wj is clearly inferior, the random 
ordering of items appears to become less and less favorable 
in comparison to the decreasing order. 

The results of the largest instance, shown in Table IV, 
confirm that the ordering of the items becomes more influential 
with increasing size of the instance. The best results have been 
obtained assigning the items in decreasing order of wj while 



TABLE II 

Identified vectors (zi , 22) for the instance with n = 200. The 

BEST OUTCOMES ARE HIGHLIGHTED BOLD. 



Item order 


Best-Fit 


Random-Fit 


Decreasing Wj 


(43, 1.000) 
(42, 1.214) 

(41, 1.927) 


(43, 1.000) 

(42, 1.881) 
(41, 1.902) 


Increasing Wj 


(50, 1.000) 
(49, 1.163) 
(48, 1.563) 


(50, 1.000) 
(49, 1.531) 


Random 


(43, 1.000) 
(42, 1.381) 
(41, 1.927) 


(43, 1.000) 
(42, 1.929) 
(41, 2.000) 


TABLE III 



Identified vectors (21,22) for the instance with n = 500. The 

BEST OUTCOMES ARE HIGHLIGHTED BOLD. 



Item order 


Best-Fit 


Random-Fit 


Decreasing Wj 


(102, 1.000) 

(101, 2.020) 


(102, 1.000) 
(101, 1.911) 


Increasing Wj 


(127, 1.000) 

(125, 1.136) 
(124, 1.177) 
(123, 1.325) 


(127, 1.000) 
(126, 1.135) 
(125, 1.152) 


Random 




(104, 1.000) 
(103, 1.951) 




(102, 1.000) 

(101, 1.972) 


(102, 1.971) 
(101, 3.168) 



the random ordering turned out to be comparably weak. When 
comparing Best-Fit and Random-Fit, Best-Fit appears to lead 
to more efficient outcomes, but still Random-Fit is able to 
identify solutions that have not been found by the Best-Fit 
algorithm. 

TABLE IV 

IDENTIFIED VECTORS (21 , 22) FOR THE INSTANCE WITH n = 1000. THE 
BEST OUTCOMES ARE HIGHLIGHTED BOLD. 



Item order 


Best-Fit 


Random-Fit 


Decreasing Wj 






(205, 1.000) 




(203, 


1.000) 






(202, 


1.287) 


(202, 1.906) 








(201, 1.910) 


Increasing Wj 


(250, 


1.000) 


(250, 1.000) 




(249, 


1.177) 


(249, 1.639) 




(248, 


1.194) 




Random 


(205. 


1 .000) 


(205, 1.000) 




(204, 


1.140) 


(204, 1.971) 




(203, 


1.897) 


(203, 2.739) 




(202, 


1.965) 


(202, 3.238)) 




(201, 


2.781) 





Common to the results of all investigated instances is that 
the approximation algorithms have not been able to identify 
the minimal solution for zi. The best solutions with respect to 
zi are still one bin larger than the minimum possible value. 

V. Conclusions 

The article presented a study on the multi-objective bin 
packing problem. We considered the objective of minimizing 
the number of bins as well as the objective of minimizing 
the average heterogeneousness of the bins, based on nominal 



attributes of the items. The two conflicting objectives led to the 
formulation of the problems as a vector optimization problem. 

A modified Best-Fit approximation algorithm has been 
presented to compute an approximation of the set of efficient 
solutions. The procedure allows the controlled consideration 
of the heterogeneousness of the bins, integrating the parameter 
in the selection process of the bins in a randomized fashion. 

Experimental investigations have been carried out on a set 
of benchmark instances, and comparison results have been 
obtained from a Random-Fit heuristic. The results are en- 
couraging, as very close approximations to the set of efficient 
outcomes have been identified. It has become clear, that with 
growing size of the instances, measured by the number of 
items n, the processing order of the items plays an increasingly 
important role. For small instances, a random order of the 
items turns out to be feasible. For large instances however, an 
order of decreasing wj is necessary to obtain good results. 

In conclusion, the presented multi-objective Best-Fit algo- 
rithm led to satisfying solutions. The running times of the 
approach remained on an Intel Pentium IV 1.8 GHz processor 
within a few seconds for each test run. We conclude that 
the algorithm may also be beneficial when computing a first, 
qualitatively good approximation of the Pareto-set which is 
then used in a more complex improvement (meta-)heuristic. 
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